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Studies have investigated the relationship between genetic variants and risk of gestational diabetes mellitus 
(GDM). However, the results remain inconclusive. The aim of this study was to investigate the association of 
rsl0830963 and rsl387153 variants in melatonin receptor IB (MTNRIB) and rsl801278 variant in insulin 
receptor substrate 1 (IRSl) with GDM susceptibility. Electronic database of PubMed, Medline, Embase, and 
CNKI (China National Knowledge Infrastructure) were searched for relevant studies between 2005 and 
2014. The odds ratio (OR) with its 95% confidence interval (CI) were employed to estimate the association. 
Total ten case-control studies, including 3428 GDM cases and 4637 healthy controls, met the inclusion 
criteria. Our results showed a significant association between the three genetic variants and GDM risk, 
rsl0830963 with a P-value less than 0.0001, rsl387153 with a P-value of 0.0002, and rsl801278 with a 
P-value of 0.001. Furthermore, all the genetic models in these three polymorphisms were associated with 
increased risks of GDM as well (P< =0.009). In conclusion, our study found that the genetic polymorphisms 
rsl0830963 and rsl387153 in MTNRIB and rsl801278 in IRSl were associated with an increased risk of 
developing GDM. However, further studies with gene-gene and gene- environmental interactions should be 
considered. 



Gestational diabetes mellitus (GDM), defined as glucose intolerance with onset or first recognition during 
pregnancy, is one of the most common medical problems and a growing public health concern^'^. It 
causes by an increase in the insulin resistance, and the condition is also aggravated by insulin secretion 
from the P cells of the pancreas^. GDM affects 1-14% of all pregnant women depending on the population 
studied^: it complicates about 1-3% of all pregnancies in the western world^, whereas 5-10% among Asian 
women^. GDM is often more common in populations with a high frequency of type 2 diabetes (T2D)^. 
Furthermore, it increases risk of adverse pregnancy outcomes and has substantial long-term adverse health 
impacts on both mothers and their offspring. Though the WHO current guidelines for GDM were published 
in 1999 and are widely used worldwide^, to date, there is still no universal recommendation for the ideal approach 
for screening and diagnosis of GDM. Thus, identif)^ing patients at a higher risk of GDM has become an important 
goal. 

Recently, extraordinary progress was made in identifying susceptible genes of complicated diseases through 
genome- wide association strategy^'^". Melatonin receptor IB (MTNRIB) and insulin receptor substrate 1 (IRSl) 
were two of diabetogenic genes associated with the developing of GDM. Melatonin is a circulating hormone 
secreted mainly from the pineal gland^^ and acts mostly through G-protein-coupled plasma membrane recep- 
tors^^. MTNRIB, located on human chromosome llq21-22^^, is a member of the G-protein-coupled receptor 
family, and one of the functional and high- affinity melatonin membrane receptors MTNRIB is expressed in 
human and rodent pancreatic islets. Studies have shown that MTNRIB is a novel candidate gene for T2D^^. 
Variants rsl0830963 (C/G) and rsl387153 (C/T) in MTNRIB have been shown with an increased risk of 
developing T2D^^. They may have a possible link in the etiology and pathophysiology of GDM. IRSl gene. 
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1 07 references initially identified 



58 abstracts screened 



49 records excluded after title review 



^ 38 records excluded at abstract review 



20 full-text articles comprehensively 
assessed against inclusion criteria 



10 references included in the meta-analysis 
9 in English/1 in Chinese 



Excluded: 

Data not extractable (n=3) 
Without useful data (n=7) 



Figure 1 | Flow chart of literature screening. 



located on chromosome 2q36^^, is expressed in insulin-sensitive tis- 
sues. It is an endogenous substrate of the insulin receptor^^, and plays 
a crucial role in the insulin signaling pathway. The IRSl gene variant 
rsl801278, a nucleotide T/C substitution in codon 972 (Gly972Arg), 
has been identified to be associated with increased risk of T2D and 

Although numerous studies have demonstrated the association 
between genetic polymorphisms and the developing of GDM, incon- 
sistent results were presented for each polymorphism among study 
populations. The purpose of this meta-analysis is to summarize the 
existing evidence on the prevalence of the genetic polymorphisms in 
patients diagnosed with GDM. 

Results 

Study selection and characteristics. The electronic database search 
identified 107 references. After applying the inclusion criteria, 10 
articles including 3428 GDM cases and 4637 healthy controls were 



ultimately included in the systematic review and meta-analysis. The 
study selection process is shown in Figure 1. 

All the 10 reports, one in Chinese^° and nine in English^^"^^, 
included cases and controls from 7 countries concerning 3 genetic 
variants in 2 genes (MTNRIB and IRSl). The detailed characteristics 
of the studies included were shown in Table 1. The distributions of 
genotypes and alleles in the individual studies were presented in 
Table 2. 

Association between MTNRIB rsl0830963 variant and GDM. For 

MTNRIB rsl0830963 variant, five studies, containing 2122 GDM 
cases and 2664 healthy controls, were included. The results of each 
allele and genetic models in this meta-analysis were listed in Table SI . 
The heterogeneity between studies was assessed, and the fixed- effects 
model and the random-effects model were employed for calculating 
the pooled odds ratio (OR). Overall, this meta-analysis showed that 
the frequency of MTNRIB rs 10830963 G allele is higher in GDM 
patients than that in the healthy controls (48.4% vs. 42.3%), and 



Table 1 | Main characteristics of studies included in this nneta-analysis 
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OGTT, the oral glucose tolerance test; ADA, the American Diabetes Association; lADPSG, the International Association of Diabetes in Pregnancy Study Groups; 

EASD-DPSG, the European Association for the Study of Diabetes-Diabetic Pregnancy Study Group; IWCGDM, International Workshop-Conferences on Gestational Diabetes Mellitus. 
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demonstrated a statistically significant positive association between 
the risk factor G allele carriers and GDM susceptibility [OR=1.24, 
95% confidence interval (CI) = 1. 14-1.35, P<0.00001)], as shown in 
Figure 2. This significant association was found in other genetic 
models as well in a fixed-effects model (GO vs. CC: OR=1.53, 95% 
CI= 1.30-1.80, P<0.00001; GG+GC vs. CC: OR=1.30, 95% 
CI= 1.14-1.47, P<0.0001; GG vs. GC+CC: OR=1.37, 95% 
CI = 1.19-1.57, P<0.0001). As shown in Figure 3. 

Association between MTNRIB rsl387153 variant and GDM. Two 

studies including 986 cases and 1070 controls focused on the 
relationship between rs 1387 153 variant and GDM. The frequency 
of the T allele was higher in GDM cases than that in controls (48.9% 
vs. 43.0%). As shown in Figure 4, our result demonstrated that the T 
allele had a positive relationship between rs 1387 153 variant and 
GDM risk (OR=1.26, 95% CI= 1.12-1.43, P = 0.0002). Various 
genetic models also demonstrated that the T allele was associated 
with an increased risk of GDM (TT vs. CC: OR= 1.56, 95% CI = 1.23- 
1.99, P=0.0003; TT+TC vs. CC: OR=1.33, 95% CI= 1.10-1.61, 
P=0.003; TT vs. TC+CC: OR=1.36, 95% CI=1. 11-1.68, 
P= 0.003) (Figure 5). No significant heterogeneity was found 
between these two studies (P=0%). 

Association between IRSl rsl801278 variant and GDM. The 

association between rs 180 1278 and GDM has been examined in 
five studies, including 1306 GDM cases and 1973 controls. Our 
meta-analysis of these studies showed that the frequency of the T 
allele of rsl801278 was higher in GDM than that in controls (8.7% vs. 
5.1%), and indicated a significant association with an increased risk 
of GDM (OR=1.42, 95% CI = 1. 15-1.75, P = 0.001) (Figure 6). As 
shown in Figure 7 and Figure SI, the dominant model and 
recessive model were also significant with GDM susceptibility. 



respectively (TT+TC vs. CC: OR=1.54, 95% CI = 1.02-2.32, 
P=0.04; TT vs. TC+CC: OR=3.01, 95% CI= 1.38-6.56, P = 0.006). 

Sensitivity analysis and publication bias. The influence of each 
study on the overall meta-analysis estimate was assessed by 
eliminating one study at a time, respectively. The OR was not 
significantly influenced by omitting any single study. 

Begger's funnel plot was used to identify individual studies in 
relation to their respective standard deviation, as shown in 
Figure 8 and Figure S2, which revealed no evidence of asymmetry. 
Egger's test was employed to provide further statistical evidence, 
similarly, no significant publication bias was found for all these three 
polymorphisms (P=0.263 for MTNRIB rsl0830963, P=0.378 for 
MTNRIB rsl387153, P = 0.149 for IRSl rsl801278). Thus, there does 
not appear to be a publication bias risk in the meta-analysis. 

Discussion 

GDM is usually recognized as a temporary form of diabetes that 
occurs during pregnancy, and is associated with an increased risk 
of complications during pregnancy and birth^°. Women with GDM 
are at a high risk of developing T2D later in life^\ and the risk of 
developing type 1 diabetes (TID) is also increased^^. Moreover, GDM 
increases the risk of macrosomia and caesarean delivery^^. Therefore, 
there is an urgent need to study the pathogenesis and establish dia- 
gnosis criteria for GDM. 

Many studies have shown that gene polymorphisms could provide 
insight into underlying pathogenetic mechanisms and the relation- 
ship between candidate genes and complex diseases. Functional 
studies showed that those diabetogenic genes took part in many steps 
of the process of developing GDM. For instance, impaired p-cell 
function (MTNRIB), insulin resistance (IRSl), and abnormal util- 
ization of glucose. In our meta-analysis, we found that the frequency 
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Figure 2 | Forest plot on the association for allelic model (G vs. C) of MTNRIB rs 10830963 and risk of GDM in a fixed-effects model. 
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Figure 3 | Forest plot on the association for the dominant model (GG+GC vs. CC) of MTNRIB rsl0830963 and GDM in a fixed-effects model. 
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Figure 4 | Forest plot on the association for allelic model (T vs. C) of MTNRIB rs 1387 153 and GDM risk in a fixed-effects model. 



of G allele in rsl0830963, and T alleles in both rsl387153 and 
rs 180 1278 respectively, were higher in GDM cases than that in 
healthy controls, demonstrating strong statistical association with 
an increased risk of GDM. 

GDM is associated with both insulin resistance and an impaired 
insulin secretion^^. GDM could develop when a genetic predisposi- 
tion of pancreatic islet P-cell impairment is unmasked by an 
increased insulin resistance during pregnancy^^. MTNRIB variants 
are related to insulin secretion and impaired p-cell function. Liao et 
al. have showed that MTNRIB is likely to be involved in the regu- 
lation of glucose homeostasis during pregnancy^^. MTNRIB 
rs 10830963 has been shown to influence the fasting plasma glucose 
(FPG)'" and to be associated with T2D''; rsl387153 has been 
reported to be associated with an increased FPG and a higher risk 
of T2D'^ In our study, we found that the SNPs rsl387153 and 
rs 10830963 in MTNRIB occur more frequently in women with 
GDM than in normal pregnant women, supporting a potential asso- 
ciation of these polymorphisms with an increased risk of developing 
GDM. This may due to the observation that MTNRIB down-regu- 
lates GCK expression and glucose- stimulated insulin secretion by 
lowering intracellular cAMP level^^. An increased expression of 
MTNRIB on P- cells leads to impaired insulin secretion. Previous 
studies have shown that the G allele of rs 10830963 polymorphism 
in the MTNRIB exhibits a higher expression of this melatonin recep- 
tor on the p-cell as compared with that of the C allele^\ 



IRSl, a substrate of the insulin receptor tyrosine kinase and a 
participant in insulin signaling^^, is related to insulin resistance. It 
plays a crucial role in the signal transduction pathway^^. 
Epidemiological studies confirmed that the prevalence of GDM is 
in direct proportion to the prevalence of T2D. A meta-analysis con- 
ducted by Jellema et al. has shown that carriers of the R972 variant of 
the IRSl gene are at a 25% increased risk of having T2D compared 
with non-carriers^^. While Morini et al. investigating 32 studies 
found that the relatively infrequent R972 variant was not signifi- 
cantly associated with T2D^^ Our result showed a significant asso- 
ciation between IRSl rs 180 1278 polymorphism and GDM risk. IRSl 
protein is expressed in many insulin -sensitive tissues, and its tyrosine 
phosphorylation can elicit the downstream effects of insulin, such as 
activation of phosphatidylinositol 3-kinase (PI3K) and translocation 
of glucose transporter 4^^. Previous studies have shown that the IRSl 
G972R polymorphism, which reduces insulin content and impairs 
insulin secretion in isolated human islets, is associated with impaired 
P-cell function^^. Evidence suggests that susceptibility to GDM has a 
genetic component, family studies indicate that GDM aggregates 
within families and is associated with a history of T2D^^. 
Several limitations were presented in this meta-analysis. Firstly, the 
number of studies included was relatively small. For MTNRIB 
rsl387153 variant, only two studies were included. Secondly, studies 
were mainly focused on Asian populations or Caucasian populations, 
other populations should also be included. Thirdly, these poly- 
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Figure 5 | Forest plot on the association for the dominant model (TT+TC vs. CC) of MTNRIB rsl387153 and GDM in a fixed-effects model. 
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Figure 6 | Forest plot on the association for allelic model (T vs. C) of IRSl rsl801278 and GDM risk in a fixed-effects model. 
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Figure 7 | Forest plot on the association for the dominant model (TT+TC vs. CC) of IRSl rs 180 1278 and GDM in a random- effects model. 



morphisms may interact with other risk factors which should be 
considered. Fourthly, the selected studies could be more subject to 
bias and artifact than prospective studies. 

In conclusion, our meta-analysis demonstrated that genetic poly- 
morphisms rsl0830963 and rsl387153 in MTNRIB and rsl801278 
in IRSl were associated with an increased risk of developing GDM. 
However, further studies with large sample sizes and accounting for 
the interaction of genetic and environmental risk factors are needed 



to understand associations between the genetic polymorphisms and 
risk of GDM. 



Methods 

Identification and eligibility of relevant studies. A comprehensive literature search 
was conducted for relevant articles published between January 2005 and March 2014 
using the electronic database of PubMed, Medline, Embase, Wanfang and CNKI 
(China National Knowledge Infrastructure). We retrieved the related articles using 
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Figure 8 | Funnel plot on the association for allelic model (G vs. C) of MTNRIB rsl0830963 and risk of GDM in a fixed-effects model (P =0.263 for 
Egger's test). 
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the following terms: "gestational diabetes mellitus", "melatonin receptor IB or 
MTNR IB", "insulin receptor substrate 1 or IRSl", "polymorphisms or variants" as 
well as their combinations. The corresponding Chinese terms were used in the 
Chinese library. References of retrieved articles were searched with no language 
restrictions. The search was focused on studies that had been conducted in humans. 
Only full-text articles and the most recent studies were included in this meta-analysis. 

Criteria for inclusion. The inclusion criteria were as follows: 1) the paper should be 
case-control or cohort studies; 2) identification of gestational diabetes mellitus cases 
was confirmed pathologically and the controls should be non-diabetic; 3) each study 
included at least one of the three polymorphisms, rsl0830963 and rsl387153 in 
MTNRIB, rs 180 1278 in IRSl; 4) genotype distribution information and OR with its 
95% CI were available; and 5) genotype distribution of control for a certain 
polymorphism must be in Hardy- Weinberg equilibrium. 

Data extraction. Two investigators independently assessed the quality of the 
included studies according to the descriptions provided by the authors of the included 
studies. Any disagreement was subsequently resolved by discussion with a third 
author. The following information was extracted from each article: first author, year 
of publication, country, ethnicity, mean age, body mass index (BMI), total numbers, 
definition and genotype distributions in GDM cases and controls. 

Statistical analysis. The overall association between genetic polymorphisms and 
GDM risk was measured by OR and its 95% CI. The Z test was employed to determine 
the significance of the pooled ORs, and a P value less than 0.05 was considered 
statistically significant. For rsl0830963, the allelic model (G vs. C) and genotype 
genetic models (co-dominant effects: GG vs. CC; dominant effect: GG+GC vs. CC; 
and recessive effect: GG vs. GC+CC) were examined; for rsl387153 and rsl801278, 
the allelic model (T vs. C) and genotype genetic models (co-dominant effects: TT vs. 
CC; dominant effect: TT+TC vs. CC; and recessive effect: TT vs. TC+CC) was 
identified. The P test was used to assess the proportion of statistical heterogeneity and 
the Q-statistic test was used to define the degree of heterogeneity. A P-value less than 

0. 10 for the Q-test and P more than 50% was considered significant among the 
studies. Data were combined using both a fixed-effects model (the inverse variance- 
weighted method) and a random -effects model (DerSimonian and Laird method)***'^^. 
The fixed-effects model is used when the effects are assumed to be homogenous, while 
the random- effects model is used when they are heterogenous. The evidence of 
publication bias was assessed by visual funnel plot inspection. Egger's regression test 
was also conducted to identify study effects (P-value less than 0.10 was considered 
significant). To evaluate whether our results were influenced by the presence of any 
individual study, we conducted a sensitivity analysis by systematically removing each 
study and reassessing the significance of the result. Statistical analyses were conducted 
in Review Manager (RevMan version 5.2, the Cochrane Collaboration, Oxford, 
England; available at: http://ims.cochrane.org/revman). All the tests were two-sided. 
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